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1 Introduction 

Gauge theories provide the best known description of the fundamental forces 
in nature. At very short distances however, physics is not known, and it seems 
unlikely that spacetime is a perfect continuum down to arbitrarily small scales. 
Indeed, physicists have started to learn in recent years how to formulate field 
theory on quantized, or noncommutative spaces. 

However, most attempts to (second-) quantize these field theories using 
modifications of the conventional, perturbative methods have failed up to now. 
It seems therefore worthwhile to try to develop new techniques for their quan- 
tization, trying to take advantage of the peculiarities of the non-commutative 
case. There is indeed one striking feature of some non-commutative gauge 
theories: they can be formulated as (multi-) matrix models. 

I will explain here the main ideas of pQ, where such a matrix formulation 
of (pure) U(n) gauge theory on the fuzzy sphere has been used to calculate 
its partition function in the commutative limit. This is done using matrix 
techniques which cannot be applied in the commutative case. 



2 The model 

Consider the matrix model for 3 hermitian N x N matrices Bi, with action 

2 / N 2 — 1 \ 
S(B) = —Trl^B^ —f + (B % + ie t3k B3B k )(B l + ie lrs B r B s )) (1) 

where g is the coupling constant and N is a (large) integer. This action de- 
scribes pure gauge theory on the fuzzy sphere S%, cp. 0E1- It is invariant 
under the U(N) gauge symmetry acting as 

Bi -> U^BiU. 
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To see that this corresponds to the usual U(l) local gauge symmetry in the 
classical limit, we first note that the absolute minimum (the "vacuum") of the 
action is given by 

Bi = Xi = 7Tjv(Ji) 



up to gauge transformation, where ttn is the ./V-dimensional representation 
of su(2) with generators Ji. Upon rescaling Xi — Xi \/ , on finds the 
generators Xi of the fuzzy sphere j3j which satisfy 



\- _ . r , _ . / 4 

This means that the vacuum of this matrix model is the fuzzy sphere. We can 
now write any field ( "covariant coordinate" ) as 

B i = X i + A i . (2) 

Then 

B l + ie ikl B k Bi = ^e lkl F kh F kl := i[X k , At] - i[X h A k ] + i[A k , A] + e Um A m . 

Notice that the kinetic terms in the field strength F k i arise automatically due 
to the shift The U(N) gauge symmetry acts on A4 as 

Ai -> U~ l AiU + U- 1 [X i ,U] (3) 

which for U — exp(ih(x)) and N — > 00 becomes the usual (abelian!) gauge 
transformation for a gauge field. One can furthermore show that the "radial" 
field 

cp := X l A t 

decouples in the large N limit, and jjTr — > J. Hence the model reduces to 
the usual U(l) Yang-Mills action 

S = ^2 J F mnF mn (4) 

in the commutative (=large N) limit. 
Monopoles. 

One can find explicitly new, non-trivial solutions of this model using the ansatz 

HI 

B t = a m x\ M) (5) 

for suitable normalization constant a m . Here x\ M ^ = itM{Ji) is the generator 
of the M-dimensional irrcp of su(2), which can be embedded in the config- 
uration space of N x N matrices if m = N — M > 0. It turns out that (J5J 



Gauge Theory on the Fuzzy Sphere and Random Matrices 



3 



describes monopole solutions with monopole charge to, and the corresponding 
gauge potential can be calculated explicitly P . Notice that negative monopole 
charges m < can also be obtained by admitting matrices B t of size TV' x TV' 
with TV < TV' <C 27V, while keeping the action as it is. 

2 

Finally, it should be noted that the correct action S — jpj of the above 
monopoles is recovered only upon a slight modification of the action Q as 
indicated in pQ, which does not affect the classical limit. The reason is that the 
"empty" blocks in (J5J if embedded in TV x TV matrices give a large contribution 
due to the first term in This is certainly unphysical (it could be interpreted 
as action of a Dirac string), and can be avoided by the slightly modified action 
(78) in pp. Then the energy of all monopoles is correctly reproduced in the 
commutative limit TV — > oo. All this extends immediately to the non-abelian 
case: 

Non-abelian case 

This model is readily extended to the nonabelian case by using matrices of 
size nTV, i.e. B t = B it0l t a = A; t° + A i;0 t° + A i>a t a where t a denote the Gell- 
Mann matrices of su(n). The action then reduces to the usual U(n) Yang-Mills 
action 

S= 7 l {Fmn '° Fmn '° + F ^n,aF mn ' a ) (6) 

in the commutative limit. Again, all "instanton" sectors are recovered if one 
admits matrices of arbitrary size M m nTV for the above action. 

3 Quantization 

The quantization of U(n) Yang-Mills gauge theory on the usual 2-sphere is 
well known, see e.g. In particular, the partition function and correlation 

functions of Wilson loops have been calculated. Our goal is to calculate the 
partition function for the YM action (JJJ on the fuzzy sphere, taking advantage 
of the formulation as matrix model. This can be achieved by collecting the 3 
matrices Bi into a single 2M x 2M matrix 

C = Co + (7) 

The main observation is that the above action can be rewritten simply as 

S(B) = TrV(C) 

imposing the constraint Cq = |, for the potential 

m = ^(C 2 -(y) 2 ) 2 

Then we proceed as 
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Z = J ' dBi exp(-S(B))) 



dC S(C - \) exp(-TrV(C)) 



dA i A i {A i ) exp(-TrV(A)) J dUS((U- L AU) - -) 

where dll is the integral over 2M x 2M unitary matrices, C = U~ 1 AU, and 
A(Ai) is the Vandermonde-determinant of the eigenvalues Ai. Here 5 (Co — |) 
is a product over M 2 delta functions, which can be calculated by introducing 

J= Co k) =k where K 

is a N x N matrix. Then 

6{(IJ- X CU) - \) = J dKexpiiTriU-^C - \)UJ)). 
By gauge invariance, the r.h.s. depends only on the eigenvalues A4 of C. Hence 
Z = J dK J dA, L A 2 {Ai)exv{-TrV{A)) J dU ex.p{iTr{U~ l AUJ -^J)) 

dK Z[J] e-i TrJ (8) 



where 

Z[J] := J dC exp(-TrV(C)+iTr{CJ)) (9) 

depends only on the eigenvalues Ji of J. Diagonalizing K = V~ x kV ', we get 

Z = J dkiA 2 (k) J dAiA 2 (Ai)exp{-TrV(A)) / dU exp{iTr{U~ 1 {A-]^)U J)) 

where J dV was absorbed in J dU. The integral over J dU can now be done 
using the Itzykson-Zuber-Harish-Chandra formula [S] 



J dU exp^TriU^CUJ)) = const 



det{e lAiJ i) 



' A{Ai)A{JjY 

which also depends only on the eigenvalues of J and C. 

In this step the number of integrals is reduced from TV 2 to 2N. This basi- 
cally means that the integral over fields on Sfj is reduced to the integral over 
functions in one variable. This is a huge step, just like in the usual matrix 
models. The constraint however forces us to evaluate in addition the integral 
over hi, which is quite complicated due to the rapid oscillations in det(e iAiJj ); 
note that Ai w ±y- Nevertheless, the integrals can be evaluated for large N 
PP, with the result 

/oo 2 
d Kl ...d Kn A 2 {k) e lK * m * exp(-|-5> 2 ). 

If 1,1 -f. . .-\-ItL n — '111 ~°° 
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Here we consider matrices of size M = nN — m, which corresponds to the 
monopole sectors with total U(l) charge m = mj + ... + m n . This can be 
rewritten in the "localized" form as a weighted sum of saddle-point contribu- 
tions, as advocated by Witten 

Z m = P(mi,9) exp(-^2^m l 2 ) 

mi + ...+m n — 771 

where P(rrii,g) is a totally symmetric polynomial in the rrii which can be 
given explicitly. In order to include all monopole configurations, we should 
simply sum over matrices of different sizes M = nN — m, for the same action 
given by V(C). One can indeed find corresponding saddle-points of the action 
which have the form 



Ai = 



( miAt 
m 2 Ai 



V o 



\ 





where 



Aj = r x A. ! 



m I 




2 1 + x 3 

becomes the usual monopole field for large N, and action becomes 



S(C< 



mi ,. . . ,m n ) \ 



2.9 



n 2 ' 



for large N. This is a standard result on the classical sphere. 

Hence the full partition function is obtained by summing 1 over all Z m , 

Z = J2z m = £ I d Ki A 2 (k) exp(-|- £ k 2 ). 



mi ,. . .,m„=-co 



Using a Poisson resummation, this can be rewritten in the form 
Z= £ AHp) exp(-27rVErf) 



or equivalently 



Z = J2 {d R f exp(-47r 2 3 2 C 2i? ). 



1 the relative weights of Z m for different m is not determined here. However, it 
could be calculated in principle 
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Here the sum is over all representations of U{n), cLr is the dimension of the 
representation and C2R the quadratic casimir. This form was found in [2] for 
the partition function of a U (n) Yang-Mills theory on the ordinary 2-sphere. 

We see that the limit N — > 00 of the partition function for U (n) YM on the 
fuzzy sphere is well-defined, and reproduces the result for YM on the classical 
sphere. This strongly suggests that the same holds for the full YM theory 
on the fuzzy sphere, and that there is nothing like UV/IR mixing for pure 
gauge theory on Sfj. This is unlike the case of a scalar field, which exhibits a 
"non-commutative anomaly" |7] related to UV/IR mixing. 
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